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Abstract 

This paper generalizes earlier work on Hamiltonian boundary terms by omitting the 
requirement that the spacelike hypersurfaces S t intersect the timelike boundary B 
orthogonally The expressions for the action and Hamiltonian are calculated and the 
required subtraction of a background contribution is discussed. The new features of 
a Hamiltonian formulation with non-orthogonal boundaries are then illustrated in 
two examples. 
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1 Introduction 



There has recently been a renewed interest, |TJ, 0, in the Hamiltonian formulation 
of general relativity, primarily motivated by the desire to extend its range of ap- 
plicability and utility to include more general boundary conditions, such as those 
which arise in black hole pair creation ||. To this date, however, it has been gener- 
ally assumed, given the standard 3+1 decomposition of spacetime, that the normal 
to the spacelike hypersurfaces is orthogonal to the normal of the boundary of the 
spacetime (although the effect of non-orthogonal boundaries has been considered, 
for example in where the action is modified in order to account for joints or 

corners where the boundaries may be nonorthogonal). The purpose of this paper 
is to derive the gravitational action and Hamiltonian, including all the terms which 
arise from the non-orthogonality of the boundaries, and illustrate the utility of such 
a derivation by calculating two examples. 

Let (M,g) be a sufficiently well-behaved four- dimensional spacetime, admitting 
a scalar time function t(x ti ), from M. onto [0, 1], which foliates Ai into a family of 
spacelike hypersurfaces, {S t }, of constant t. The boundary of M. consists of the 
initial and final spacelike hypersurfaces So and Si, as well as a timelike boundary 
B, hereafter called the three-boundary For each S 4 , we can define a two-surface, 
Bf — St n B, which bounds the hypersurface. The family of two-surfaces {B t } 
then foliates the three-boundary B. The spacetime and its submanifolds are shown 
in figure 1. The tensors defined on the surfaces are given in table 1, and are an 
amalgamation of the notation adopted in ]I[ and 0. Greek letters are used for 
indices on A4, while middle roman letters (i ■ ■ -p) are used for indices on S f , middle 
roman letters with a circumflex for indices on B, and early roman letters (a ■ ■ ■ d) are 
used for indices on B t . Tensors on any of the submanifolds can also be considered as 
tensors in A4, by the obvious embedding, and in this context are denoted by greek 
indices. However, care must be taken with raising and lowering indices. An index is 
considered to be lowered or raised by the metric corresponding to the type of index 
used. 

We will work in a coordinate system adapted to the time function and the three- 
boundary, that is, the first coordinate is t, and B is a surface of constant x l . This 
allows us to write the metric g^ u in the usual ADM |J decomposition, 

ds 2 = -N 2 dt 2 + h lj (dx i + N i dt)(dx j + N j dt), (1) 

where iV is the lapse function (taken to be the positive square root), and N l , which 
lies in the tangent plane of St, is the shift vector. We can define a timelike vector field 
t^, which we interpret as connecting corresponding points on adjacent hypersurfaces, 

t n = Nn n + N n^ (2) 

where n M is the unit normal to the spacelike hypersurfaces. Thus, given a point on a 
hypersurface, its evolution normal to the hypersurface is governed by the lapse func- 
tion, while the shift vector dictates its evolution tangent to the hypersurface. The 
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Figure 1: The spacetime manifold M. and its sub manifolds are shown, with one 
spacelike dimension suppressed. The unit normals n M , and to the spacelike hy- 
persurface E t and the timelike boundary B are shown at points on the two-surface 
B t . 



orientations of the unit normals n^ 1 and are fixed by requiring that they be future 
pointing and outward pointing, respectively. The intersection of the two hypersur- 
faces is not required to be orthogonal, since this would be an unnecessary restriction 
on the spacetime. The non-orthogonality is measured by using the variable 



(3) 



which clearly vanishes in the orthogonal case. 

If we want B to be mapped into itself by i M , then clearly ft 1 must lie in the three- 
boundary, and hence it must be orthogonal to the normal u^. For this to occur, it 
is necessary that the lapse function and shift vector satisfy 
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(4) 
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Table 1: The naming conventions for the tensors on the spacetime Ai, and the 
hypersurfaces and surfaces embedded therein. 



But, more generally, it is not necessary for £ M to lie in the boundary, and indeed it 
is impossible if you are considering a spatial translation. 

In section |2] of the paper, the action and Hamiltonian are calculated for metrics of 
Lorentzian signature. In order to obtain a finite answer for a non-compact spacetime 
it is necessary to consider the action and Hamiltonian relative to a background 
spacetime jlOfl . This process is outlined in section [|. In section |, two examples are 



presented in order to illustrate some of the properties of the Hamiltonians obtained 
from spacetimes with non-orthogonal boundaries. The appendices contain various 
kinematical equations, some relations based on the particular coordinate system 
adopted in this paper, and a summary of the corresponding Euclidean results. 

2 Calculation of the Action and Hamiltonian 

The standard Hilbert action for Lorentzian general relativity is 

I[M,g] = — I d 4 x v ^gR+-r i d 3 xVhK+-[d 3 x v ^Q 
2k Jm k Jt, k Jb 

1 r B i 

+ — / d 2 x ^/o-sinh^ 1 ^, (5) 

K JB 

where k is 8nG, the integral between E and Ex is notation for the integral over 
the final hypersurface, E 1; minus the integral over the initial hypersurface, E , and 
similarly for the integral over the initial and final two-surfaces, Bq and B\, (which 
bound the initial and final hypersurfaces So and Ei). The final term, referred to 
as the corner term, is necessary in order to ensure that the variation of the action 
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with respect to the intersection angle rj vanishes, |J, and is the only effect of non- 
orthogonal boundaries which has been considered previously. In order to obtain 
the Hamiltonian, we need to factor the integrals into an integral over time t, and 
integrals over spacelike surfaces. We would also like to separate out the terms in the 
volume integral which are pure divergence, and convert them to boundary integrals. 
Let Im, Ib, It,, and Ib be the integrals over the corresponding manifolds in equation 
(§• ' 

Using equation fl83"l) , we can substitute, for the curvature scalar R, terms which 
lie in the hypersurface S 4 . Thus, the volume integral becomes 

J M = ^-[ d 4 x^[Jl + K llu K^-K 2 + 2V ll (n»K-a»)}, (6) 

2k J M 

where a M = n u 'V u n tJ- is the acceleration of the unit normal n^. We can convert the 
final term in the volume integral to a surface integral over the boundary of Ai. 
Thus, Im simplifies to 

I M = 7rl d'x^n + K^-K 2 ]. (7) 
ZK j m. 

The integral over the three-boundary becomes 

I' B = - f d 3 x y/=j[Q + uJn^K - O], (8) 

K JB 



while the integral over the initial and final hypersurfaces becomes 

1 pi 



4 = I ^ d 3 x Vh[K + n^rfK - a")]. (9) 

K J Sn 



But, since n M a M vanishes, and the K terms cancel, we see that 1^ disappears. Hence, 
the action reduces to I' M + I' B + Ib, 

I = — I dt ( d 3 x NVh[JZ + K^K^ - K 2 ] + - f d 3 x y/=i[S + rjK - u^] 

r-Bi 

+ / d 3 x ■y/as&rih~\ (10) 

J B 

We now want to maneuver the action into canonical form. Using equation (jS8|) , we 
can write I' M in terms of canonical variables, 

I' M = JdtJ^ d 3 x \P^\ IV — NH — 2P'»D IM N U ], (11) 

where 7i is the energy constraint, which vanishes on any solution of the field equa- 
tions. From equations ( |103| ) and ( 1 14jj we can write I B as 



I' B = - [dt [ d 2 xNy^[k - AVV^], (12) 
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where i> M = A(n M — r]u^) is the normalized projection of n M onto the three-boundary, 
and A = (1 + rf)~^ is the normalization constant for the unit vectors and v^. 
We can write the nonorthogonal part of the boundary integral as the sum of a total 
derivative and a second term, 

I' B =- f dt f d 2 x [N^k - v ^(V^(^sinh~ 1 ^) + sinlT^VX)]. ( 13 ) 

K J J Bt 

If we convert the total derivative to a boundary integral over Bq and B\, then it 
will cancel Ib, and the action will simply be the sum of the integrals over Ai and 
B, where the three-boundary integral is now 

I" = -[dt[ d 2 xNy^(k + Asinh^VX)- ( 14 ) 

K J JB t 

Finally, to obtain the action in canonical form, we need to remove from I' M the 
term involving the derivative of the shift vector. Using equation (|9l]), we obtain 

I' M = [dt [ d 3 x [P^'h^ -NH- Nm^ - 2D^N V % (15) 

where H.^ is the momentum constraint, which also vanishes on any solution of the 
field equations. Converting the final term into a surface integral over B t , the volume 
contribution reduces to 

I" M = [dt ( d 3 x [P^V -NH- N^W], (16) 

while the boundary term becomes 

if = - ( dt f d 3 xN^[k + Asinh-^VX];- 2 I dt I d 2 xr^ u N u , (17) 

where Pg v is the tensor density P^ v with the correct area element for B t , as given by 
equation (|104|) . Thus, by factoring out the time integral, we can express the action 
in canonical form, 

/ = Jdty^d 3 xP^h llu -(H c + H k + H t + H m )Y (18) 

where the Hamiltonian, 

H[M,g] = H c + H k + H t + H m , (19) 
is a sum of four distinct terms: 
1. a constraint term, 



f d 3 x[NH + N^Wl (20) 
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2. a curvature term, 

H k = — [ d 2 xNy/^k, (21) 

K JB t 

3. a tilting term, 

H t = —f (fxNX^sinh^rjVaV^, (22) 

K J B t 

4. and a momentum term, 

H m = 2 I d 2 xr^N u . (23) 

In the case of orthogonal boundaries, the curvature term, H^, is usually taken to 
give the mass of the system, while the momentum term, H m , gives the linear or an- 
gular momentum. However, it will be shown in section f| that the two contributions 
get mixed up in the non-orthogonal case. 

We also see that the only place in which the non-orthogonality appears explicitly 
is in the tilting term, H t . One can see that the action must depend on the way one 
chooses the angle between E 4 and the three-boundary to vary with time. If the 
resulting two-surfaces, B t , are independent of time, as they would be on an inner 
horizon, then the tilting term vanishes, since the total non-orthogonal contribution 
to the integral over the three-boundary is simply the negative time derivative of 
the corner term, and hence the two terms will cancel. If the two-surface is time- 
dependent, as it can be on the boundary at infinity, then H t will in general be 
nonzero. However, this term will cancel when we subtract the action of the back- 
ground, as detailed in the next section. Since the action must be independent of 
the way in which we choose our boundaries to intersect, we see that a Hamiltonian 
treatment of the action will not work unless we include a background subtraction 
to remove the tilting term. 

3 Background Spacetime 

If the time surfaces S t are non-compact, the action is calculated by evaluating the 
action on a compact region, and then letting the boundary tend to infinity. This is 
problematic, since the Hamiltonian will generally diverge as the boundary is taken 
to infinity. However, it makes sense to define the physical Hamiltonian, H p , to be 
the difference between the Hamiltonian for the space under consideration, and the 
Hamiltonian for some background solution of the field equations, which can be re- 
garded as a ground state for solutions with that asymptotic behaviour. Quantities 
defined on the background spacetime are indicated by a tilde. A minimum require- 
ment for a solution to be regarded as a ground state would seem to be that it had a 
timelike Killing vector, but one might also ask that it was completely homogeneous 
and had three linearly independent Killing vectors as well. The existence of the 
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timelike Killing vector is necessary for a suitable definition of energy. The usual 
background is Minkowski space, but one can consider other backgrounds, such as 



anti-de Sitter space [JTTJ] , the Robinson-Bertotti metric |]12| , or the Melvin universe 
|Hifl . The last is not homogeneous, but has been used as the background metric for 
the pair creation of charged black holes || . Choosing an appropriate background for 



Kaluza-Klein spacetimes presents additional difficulties |b|, [L5[], due to the presence 
of the compactified dimensions. These spacetimes will not be considered here, but 
will be addressed in forthcoming papers fl6|, 

The induced metrics on the three-boundary should agree in the two solutions, 
but in some cases they will agree only asymptotically, in the limit that the three- 
boundary goes to infinity.^ It is a delicate matter to choose the rate at which the 
metrics on the boundaries approach each other, and it will depend on the asymptotic 
behaviour of the lapse function and shift vector under consideration. We shall 
assume that the necessary fall-off conditions are satisfied, and shall therefore take 
the induced metrics on the boundary to agree. 

Since the background metric is taken to be a solution of the field equations, 
the physical Hamiltonian will only have a constraint contribution from the actual 
spacetime (and only if it is not a solution), 

E p c = f d 3 x[NH + N^W]. (24) 

J St 

Because the induced metrics of the actual spacetime and its background agree on the 
three-boundary, they will have the same volume element, NXy/a. Furthermore, we 
can always take the angle between E 4 and B to be the same for the two spacetimes, 
hence the respective values of A will also be the same. Thus, the curvature term is 
simply 

H? = --[ d 2 xN^(k-k). (25) 
n JB t v ' 

If the background slices are chosen such that the conjugate momentum vanishes, then 
the momentum contribution to the Hamiltonian will simply be from the spacetime 
under consideration. However, the constraint of matching i] in both spacetimes 
implies that we cannot always ensure that the momentum density vanishes in the 
background. Thus, 

E p m = l\ d 2 x(r,P^N u -r,P^N u ). (26) 

Since i] is the same for both spacetimes, and the vector lies entirely in B, the 
tilting terms will be equal for both spacetimes, and hence will cancel each other 



3 If the spacetime has a horizon, then it may be possible to analytically continue the solution 
through the horizon to obtain a second asymptotic region of space. In this case, we want to take 
Ti t to have an inner boundary on the horizon, rather than two asymptotic regions. Alternatively, 
we could use a star, or similar physical object, to eliminate the extra asymptotic region. 
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when subtracted. Thus, there will be no contribution to the physical Hamiltonian 
from the tilting term, 

Hf = 0. (27) 

The vanishing of the tilting term is necessary for the Hamiltonian formulation to be 
well-defined, because otherwise we have an unacceptable dependence of the action 
on the way in which the boundaries intersect. Note that on an inner boundary, such 
as a horizon, there may not be a subtraction, since the background may not have a 
horizon, but, as noted at the end of section 0, the tilting term will vanish anyway, 
since the horizon two-surfaces are time- independent. 

For a Killing vector, of the background spacetime, we can obtain a conserved 
charge on the spacelike hypersurfaces by decomposing £ M in terms of a lapse function 
and shift vector, and calculating the corresponding Hamiltonian. Thus, assuming 
that t M is asymptotically equal to the time translation Killing vector of the back- 
ground spacetime, the energy, E, which is the conserved charge associated with time 
translation, is simply the value of the physical Hamiltonian, H p . 



4 Examples 

4.1 Schwarzschild Spacetime with Flat Spacelike Slices 

We first want to consider a simple example which has a non-orthogonal intersec- 
tion of boundaries, but for which the tilting term H t vanishes, and where there is 
no non-trivial spatial linear momentum. We begin by computing the terms in the 
physical Hamiltonian, first for the Schwarzschild spacetime, and then for the back- 
ground Minkowski space. Two different background coordinate systems are used-a 
minimally matched system in which the induced metrics agree, but the values of 
rj differ, and a correctly matched system in which both the induced metric and r\ 
agree. Once the terms have been computed, we can obtain the physical Hamilto- 
nian. The calculations of this example and the next one both made extensive use of 
the GRTensorll package for Maple |jTB| . 



If we consider the Schwarzschild solution in standard static coordinates, 

2M. , 2 1 
— )dt 2 + - 
r 1 



ds 2 = -(1 - — )dt 2 + 2M dr 2 + r 2 dn, (21 



and then define a new time coordinate t' by 



dt = dt' — — jwdr, (29) 



then the Schwarzschild line element becomes 




ds 2 = -dt' 2 + (dr + J—dt') 2 + r 2 dtt (30) 
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= ^i-^L )d ^ + 2 J—drdt' + dr 2 + r 2 dn, (31) 

T V T 

which is the Painleve and Gullstrand line element, most recently investigated 



in 



20|j . This coordinate system is of interest primarily because it has flat spacelike 
slices and the intersection of a hypersurface of constant time and one of constant 
radius is non-orthogonal. The lapse function and the radial component of the shift 
vector are non-vanishing, 



iV=l, and N r = ^^. (32) 

As stated in the introduction, we want to consider spacelike hypersurfaces S t / of 
constant time t', and a three-boundary B of constant radius r. In order to calculate 
the physical Hamiltonian, we need to compute the Hamiltonian for a fixed three- 
boundary which we then let tend to infinity. Hence, take the fixed three-boundary 
to be the hypersurface of radius R, denoted B R . Since the calculation is independent 
of the choice of spacelike hypersurface, we do not need to fix E t / . We find that t^u^ 
vanishes, and hence t M lies in the three-boundary, so that B R is evolved into itself 
by £ M . The intersection of E t / and B R is not orthogonal, but is characterized by the 
variables 

/ 2M , / 2M , . 

^ = -V^2M' and ^f-^T (33) 

The two-surface B R = T, t > R B R is a sphere of constant t' and r. 

As noted above, the induced metric on £ t / is completely flat, and hence the 
induced metric on B R is simply the standard two-sphere metric. Thus, these metrics 
have identical volume elements, while the volume element of the three-boundary B R 
contains a contribution from the non-trivial value of A, 



Vh = B 2 sin 9 = y/a, and NX^ = y 1 - sin 9. (34) 

We now would like to calculate the terms which contribute to the Hamiltonian. 
The trace of the extrinsic curvature of the two-surface, B R , is 

(35) 

The two-surface B R is independent of time, so V^u M vanishes, and hence there will 
be no tilting contribution to the Hamiltonian. If we calculate the momentum tensor 
density on B R , P£ u , and contract it with the unit normal, r M , and the shift vector, 
N^, we obtain 

r^ v N v =—sin9, (36) 

K 

which is independent of the fixed radius R. 
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The natural choice for the background spacetime is Minkowski space, since it 
has the same asymptotic behaviour as the Schwarzschild solution. We first consider 
a coordinate system in which the induced metric on the three-boundary of radius R 
agrees with the Schwarzschild case, but the intersection angle, characterized by rj, 
is different. The correct background coordinate system, with matched r) values, is 
then used. The Hamiltonian is identical in both cases, but the contributions from 
each of the terms depends on the coordinate system chosen. 

Consider the static coordinate system with a scaled time coordinate, 

ds* = _(1 _ 2 -^-)dt 2 + dr 2 + rW. (37) 
R 

As in the Schwarzschild case, we consider a fixed three-boundary of radius R, B R . 
The induced metric on B R is then the same as the induced metric on the three- 
boundary B R , and hence we can equate the two boundaries. The lapse function and 
shift vector are not the same as in the Schwarzschild case, but instead 




while the shift vector vanishes. The disagreement in the lapse is due to the dis- 
agreement between the boundary intersection value r\ in the two solutions. Since 
the boundaries in the background coordinate system were chosen to be orthogonal, 
fj vanishes, and A reduces to one. 

We now want to calculate the terms contributing to the Hamiltonian. The trace 
of the extrinsic curvature of the two-boundary is 

fc=i (39) 

As in the Schwazrschild case, the two-surface B R is time- independent, and hence 
V M f M is zero. The momentum density P^ v also vanishes, so that 

r^ v N u = 0. (40) 

We can now calculate the physical Hamiltonian. Since both the Schwarzschild 
and Minkowski spacetimes are solutions of the field equations, the constraint term 
vanishes. The curvature contribution to the physical Hamiltonian can be calculated 
by integrating equations (|35f) and (|39|) over the two-surface B R (taking care to 



include the factors arising from the difference in A values between the two solutions) 
Hjf - lim / <!<■> I ' ' \ 1 ^^It, ' ..., -tA- Ml) 




If we expand the square root in a Taylor series, then the infinite contributions cancel, 
and we are left with only a finite value, 

H R = -M. (42) 
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The tilting contribution is zero because V M t> M vanishes in both the Schwarzschild and 
Minkowski cases. The momentum contribution is due entirely to the Schwarzschild 
term, equation (j36j). 

r 2ir r l 2M 

H? = 2 d<f> d(cos6) = 2M. (43) 

JO J -l K 

Thus, the Hamiltonian is 

H p = H[ + H? = M, (44) 

as anticipated. Note that there were contributions from both Hj? and H^, contrary 
to standard expectation. However, as will be shown, this is due to incorrect matching 
of the background. 

We would now like to consider a background coordinate system in which the 
value of r) agrees with that of the Schwarzschild solution. If we introduce a new 
time coordinate 

dt = dt' / R dr. (45) 



I __ 2M 
1 R 



then the line element (|37p becomes 



ds 2 = -(1 - -j^-W 2 + 2W-^-Wl - — dt'dr + (1 - -j^)dr 2 + r 2 dVL 2 . (46) 
The lapse function and radial shift vector are non-vanishing, 



2M 2M 

N=l, and iV- = y — yi-— , (47) 

and the lapse agrees with the lapse of the Schwarzschild case. The induced metrics 
on a three-boundary of fixed radius R still agree, but now the boundary intersection 
is non-orthogonal, since 



2M ~ / 2M 

n = -\ , and A=W1 , 48 

' V R-2M' V R 

which are equal to the Schwarzschild values. If we calculate the trace of the extrinsic 
curvature of the two-surface, Bft, we obtain 

k = . 1 (49) 

rJi ™ 



R 

The conjugate momentum is now non-trivial in this coordinate system, 

r,(P a r"N u = — sm9. (50) 

K 
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We see that just as in the other coordinate system, 'V fl v fl vanishes. 

We now want to calculate the physical Hamiltonian. To calculate the curvature 
contribution, we integrate equations (^5|) and Q49|) over the two-surface (where now 
since the values of A agree, we can use equation 



I /-27T pi 

H p = lim / d<j> d(cos9)R 2 

K R-^oo Jq 7-1 



1 



r 



(51) 



When we expand the square root in a Taylor series, we obtain 

H p = M. (52) 



Since the two momentum terms, ( p6| ) and are identical, they will cancel when 
subtracted, and hence H m will vanish. As in the previous case, V^v^ disappears 
in both the Schwarzschild and Minkowski spacetimes, and hence there will be no 
tilting term contribution to the Hamiltonian. Thus, the only contribution to the 
physical Hamiltonian is from the curvature term, 

H p = H P = M, (53) 

as anticipated. Note that this is an example where the constraint of matching rj 
between the spacetimes means that we cannot pick the background spacetime such 
that vanishes. This is contrary to the situation when only non-orthogonal 
boundaries are considered. It is obvious from the symmetry of the spacetime that 
the linear and angular momenta vanish. 

As noted in section |3], in order to avoid a second asymptotic region we must 
consider an inner boundary on the horizon, at r = 2M. There is no background 
contribution, since Minkowski space has no horizon there, and hence the result will 
be the same for both choices of background coordinates. If we calculate the curvature 
contribution at the horizon, we obtain 

H p = -2M, 

while the momentum contribution (which was noted to be independent of radius) is 

H p = 2M. 

Thus, we see that there is no net contribution from the horizon. 



4.2 Schwarzschild Spacetime with Tilted Spacelike Slices 

In this example, we consider a more complicated slicing of the Schwarzschild space- 
time, corresponding to a Lorentz time boost in the z direction, which leads to a 
non-trivial tilting term. In addition, we obtain a non- vanishing value for the spatial 
linear momentum. 
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If we take the Schwarzschild solution in static coordinates, and make the variable 
substitution, 

t' = , 1 =(t-vr cos 0) , (54) 



then the line element becomes 
2M 



ds z = - 1 - 



V 1 — v 2 dt' + w cos 6*(ir — sin 



1 _ 2M 



dr 2 +r 2 dQ. (55) 



The lapse function is given by 

AT 



r(r-2M)(l - v 2 ) 



^ r 2 _ v2(r _ 2M)(r - 2M cos 2 ) ' 
while the shift vector has nonzero components in the r and directions, 



(56) 



N r = -vVT 



2M 



cos 0, and N e = vy/T 



2M 



rsin0. (57) 



As before, we take E t / to be the hypersurface of constant t', and we fix the three- 
boundary, B R , to be the hypersurface of constant radius R. The intersection pa- 
rameter rj is 



V = 



v(r-2M) cos0 



yjr 2 - v 2 (r 2 - 2Mr{\ - cos 2 0) + AM 2 cos 2 0) 



(58) 



The two-surface 5^ is simply a two-sphere of constant radius R and time 

We now want to calculate the quantities which contribute to the physical Hamil- 
tonian. Expanding the relevant expressions in Taylor series, only terms which yield 
a non-vanishing value, in the limit as R goes to infinity, are given. To simplify the 
notation, we set x = cos 0. The extrinsic curvature of the three-boundary, Bj}, is 



k = 



{2 - v 2 [\ - x 2 ])^r^p i 



2(1 -v z ) -x 2 vH\ 



X 



M 



[l - V 2 [l - X 2 ]f/ 2 R y/T=tf(l -V 2 [l -X 2 \f/ 2 R?' 

The derivative of r\ in the direction of is 

vx(2 - [1 - m[l - x*]v>) 



(59) 



r 3 / 2 (l-[l-^][l-:r> 2 )' 



(60) 



If we calculate the conjugate momentum, and contract it with r M and the shift 
vector, we obtain 

™,„», M n v 2 (3x 2 + l) 
r„P^N u = sin#- 



M a 



2k 



(61) 



We can now define the background spacetime, and calculate its contribution 
to the physical Hamiltonian. We need to find a coordinate system such that the 
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induced metric on a boundary of constant radius R agrees with the induced metric 
due to the Schwarzschild solution. If we start with Minkowski space with scaled 
time, given by equation fl37D, and then make the coordinate transformation, 



^ttM''-WM co 4 (e2) 



where 



/■ | 2M \ 9 f, — jw 

/W = 2m(1 + a/1- — I e-V*-- (63) 



then the resulting coordinate system gives the same induced metric, value of rj, and 
lapse function on the three-boundary of constant radius R as those found in the 
Schwarzschild metric. If we calculate the trace of the extrinsic curvature of the 
two-surface, we find that 

~ _ (2 - v 2 [l - x 2 ])VT^ 1 _ (2(1 - v 2 ) - x 2 (A - v 2 [3 - x 2 ]))v 2 M 

(1 - V 2 [l - X 2 }f/ 2 R v /j-^2(l_ t; 2[ 1 _ x2 ])5/2 R2- (> > 

Thus, if we integrate the difference between the Schwarschild value, given by equa- 
tion (^9[), and this background value, then the infinite parts cancel, and we obtain 
the curvature contribution, 

H P k = -^L=. (65) 
V 1 — v z 

The conjugate momentum, contracted with the unit normal f M and the shift vector, 
is 

~ ~ M 3r 2 — 1 

2k, a/1 - v 2 

We see that if we integrate this over the two-surface it vanishes, and hence the 
momentum contribution comes entirely from integrating Schwarzschild term, given 
by equation (pi]), 

Hi ^ (67) 
V 1 — u 

Finally, we note that, by construction, the tilting terms are equal in the Schwarzschild 
and Minkowski systems, and hence cancel each other. Thus, as expected, there is 
no tilting contribution to the Hamiltonian. Hence, if we add the curvature and 
momentum values, we obtain the value of the physical Hamiltonian, 



H p = H[ + H? = MVT^ 2 . (68) 

From its static value, the Hamiltonian has been decreased by the inverse of the 
boost factor. As will be shown below, this decrease in energy is accounted for by 
a non-zero value of the linear momentum in the z direction. Note that unlike the 
case when orthogonal boundaries are used, the physical Hamiltonian now contains a 
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non-trivial contribution from the momentum term. On the horizon the shift vector, 
the tilting term and the curvature term vanish, so it provides no contribution to the 
Hamiltonian. 

We now want to consider the conserved charges arising from the Killing fields of 
the background spacetime. The asymptotic value of the background spacetime, in 
Cartesian coordinates, is 

ds 2 = -(1 - v 2 )dt' 2 - 2vVl - v 2 dt'dz + dx 2 + dy 2 + (1 - v 2 )dz 2 . (69) 

This is related to the standard static background by the transformation 

t' = —L=(t-vz). (70) 
V 1 — v 2 

We can associate four Killing fields, t^, and xf with the translation invariance of 
each of the coordinates. Since is equal to the time evolution vector the 
conserved charge associated with time translation, Vt, is simply the value of the 
physical Hamiltonian, Ma/1 — v 2 . If we consider the spatial Killing vectors xf, we 
find that by symmetry, the linear momenta in the x and y directions vanishes, but 
due to the transformation (|7D|), there is a nonzero value of the z momentum. The 
z momentum will contain only a contribution from the momentum integral. If we 
contract the conjugate momentum with the unit normal and the Killing vector, we 
obtain 

z ll P?'r v = Mv(3a? + l). (71) 
The corresponding background value is 

ZyPFfv = Mv(3x 2 - 1). (72) 

Thus, the background contribution will integrate to zero, while the Schwarzschild 
term yields 

V x = J d 2 x~z^ v r v = Mv. (73) 

The resulting energy momentum vector is 

= (MVl - v 2 , 0, 0, Mv) . (74) 

If we calculate the norm of the vector, with respect to the asymptotic background 
Minkowski space, we obtain 

= — M 2 , (75) 
indicating that it has transformed correctly as a four-vector. 
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A Kinematics 



In this appendix, various standard kinematical formulae are presented, and some 
relations used in the paper are derived. 

A.l The hypersurfaces S f 

The basic quantities which are induced on the hypersurface S t by the metric g^ u 
and the unit normal are the first and second fundamental forms, generally called 
the induced metric and the extrinsic curvature, 

K, v = h«V a n v . (77) 

Any tensor may be projected onto E t by using the projection tensor . In this 
way, we can define the induced covariant derivative on £ t as the projection (of every 
index) of the covariant derivative of the tensor in M.. For example, 

DfjT vp = h^h^h^V^. (78) 

where T vp is a tensor field defined on S 4 . We now want to write the curvature scalar 
for Ai in terms of quantities defined on E t . We begin by decomposing it in terms 
of the Einstein and Ricci tensors, 

R = 2{G lxv - R^rfn". (79) 

From the Gauss-Codazzi relations, we obtain the initial value constraint 

G^n" = \{K- K^K^ + K 2 ). (80) 
By definition, the Riemann tensor satisfies 

R^n v = -^(V^V, - VvVjn", (81) 
which, after some simplification, yields 

R^n v = K 2 — K^K^ - V^K - a"), (82) 

where a M = n v, W u n p is the acceleration of the unit normal n M . If we substitute 
equations fl5U|) and into equation ([73), we obtain the desired expression for R, 

R = TL + K^K^ - K 2 + 2V fl (n p K - a"). (83) 

We now want to write this expression in terms of the canonical variables, P^, 
h^v, N, and N^. We begin by expressing the extrinsic curvature in terms of these 
canonical variables, 

K lut = ^(V-2L> ( ^)), (84) 
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where h^ v indicates the Lie derivative along the evolution vector t^. If we substitute 
this into the action, we can calculate the momentum conjugate to the hypersurface 
metric h^, 

pw = —y/h{K ia ' - Kh? v ). (85) 
If we calculate the term due to the extrinsic curvature in equation (|83"|), 



2k kN 

K^K*" — K 2 = ——[P^h^ - -j={2P, v P^ - P 2 ) - 2P^D fl N u ], (86) 

we can then define the energy constraint as the term in the action which vanishes 
when the variation of N is set to zero, 

H = 4=(2P^ - P 2 ) ~ ^Tl, (87) 

and hence we see that 

— NVh[TZ + K^K^ - K 2 } = P^'h^ -NH- 2P^D^N U . (88) 

If we write 

P^DpNv = D^P^Nu) - N^D U P^, (89) 

then we can define the momentum constraint (which vanishes when the variation 
due to the shift vector vanishes) as 

W = -D V P^ V , (90) 

and hence we obtain 

P^D^N V = D^P^Nv) + N»Hr (91) 
A. 2 The three-boundary B 

The induced metric and extrinsic curvature of B are given by 

= 9/iu - u^u U) (92) 

6^ = ifVaUu. (93) 

Since we are not assuming that n M is in the tangent space to B, the scalar product 
of the two, 

V = ™X> ( 94 ) 
may be nonzero. Moreover, if r\ is nonvanishing, then the normalized restriction of 
n» to B, 

v » = \^ a n a = \(n» - 77O, (95) 

will not, in general, be equal to n M . If we want t> M to be a unit timelike vector, we 
find that 

A = -t=J==j. (96) 
Vl + T 
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A. 3 The family of two-surfaces B t 

We want to consider B t as a hypersurface embedded in £ t . Thus, the normal to 
B t must lie in the tangent plane to £ t . Therefore, if is not orthogonal to 
on B, then we cannot take to be the normal to B t . Instead, we must take the 
normalized restriction of to S t , 

= Xh^Ua = X(u^ + i]n^), (97) 

where the normalization constant is again A, as given by equation (Pop. 
Thus, the induced metric is given by 

<V = Kv ~ V" = 9i»> + X 2 ^^ - u^Uv - 2r]n^u v) ). (98) 

We now want to express the extrinsic curvature of B t in terms of quantities 
defined on £ t and B. By definition, it is 

= a^ a D a r u . (99) 

If we expand a^Da, and take the trace, we obtain 

k = V M r M + •nPn v V v rr (100) 

But, using the orthogonality of n M and r M , we see that 

nW v r^ = -r^V^ = -Xu v a v , (101) 

and by the definition of the extrinsic curvature of the £ t and B, 

V M r^ = A[6 + r]K + A 2 (r^ - rju^V ^7]]. (102) 

Combining the two results, and using the definition of the restriction of to B t 
from equation (p5|) , we obtain the following value for k, 

k = X[e + rjK -u^ + Xv^V^ri]. (103) 

Note that since t> M is the projection of n M onto B, the derivative in the expression 
for k is in the direction perpendicular to B t which lies in B. 

If we consider the momentum tensor density conjugate to hij, as given by equa- 
tion (Hp) , then it takes a different value when viewed as a density on B t , because 
the volume element has changed. On B t it is given by 

= —^(K^ - KbT). (104) 
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B Coordinate Conditions 

We now present some formulae due to the particular coordinate system chosen, that 
is, oriented with respect to the hypersurfaces E t , and the three-boundary B. The 
formulae are primarily aimed at calculating the volume elements in terms of each 
other. In our chosen coordinate system, we see that the unit normals are given by 

u„ = (-TV, 0,0,0), (105) 
rv = (0,^,0,0). (106) 

If we then calculate 77, we see that 

V = gTn^ = -g 01 N^. (107) 

We can use this expression to relate the lapse function and the intersection variable 
to the first component of the vector u^, 

u° = g°^ = g 01 ^=-^. (108) 

To calculate the relationships between the determinants of the metrics, we use 
the matrix identity, 

{A ~ 1)l > = dhi { - iy+ ' detA{h3) > (109) 

where A(i,j) is the matrix formed by removing the i th row and j th column from A. 
Then, since <7y(0, 0) = h^, we see that 

g 00 = ~h. (110) 

But, g 00 = —N~ 2 , and hence we obtain the familiar relation 

^g~ = N\fh. (Ill) 
Using the definition 7afe(0,0) = a a b, equation ( |109| ) yields 

7°° = (112) 

7 



If we then substitute in the value of 7 00 , 

. 1 

■JP + ~N 2 ' N 2 ^ ' '' ' N 2 X 



y» = r - «v = -(4 + = -4(1 + , 2 ) = --JL (113) 



we obtain the desired relation, 

V ^=N\^. (114) 
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C Euclidean Formulae 



There are several important definitions which are changed when we deal with a Eu- 
clidean rather than a Lorentzian metric. The method of obtaining the Hamiltonian 
from the action remains the same, but several of the final results contain negative 
signs relative to their Lorentzian counterparts. Euclidean quantities, such as the 
action and Hamiltonian, will be denoted by a circumflex. The Hilbert action for 
Euclidean general relativity is 



167T 

cB 



— ( JgR - — / Sl VhK - — [ Ji® 

6n Jm 8n Js Stt Jb 



/ d 2 x v^acos" 1 rj. (US) 
Jb 



Performing the same steps in decomposing the action as were followed in the Lorentzian 
case, we obtain 

/ = J dt { jT d 3 x P^'h^ + (H c + H k + H t + H m )Y (116) 
where the Hamiltonian, 

H[M,g] = H C + H k + H t + H m , (117) 
is again a sum of four distinct terms: 

1. a constraint term, 

H c = -[ d 3 x[NH + N^], (118) 

2. a curvature term, 

H k = --[ d 2 xN^k, (119) 

K JBt 

3. a tilting term, 

H t = --( d 2 xNXy/^ cos -1 ?|V/, (120) 

K JB t 

4. and a momentum term, 

H m = -2 I d 2 xr^N v . (121) 

JBt 

In order to derive these results, it is necessary to note that is now a spacelike 
vector, and hence it induces a metric 
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The normal to the two-boundary is then 
where the normalization constant A is now 

Apart from these changes, the analysis is identical. 
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